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1.  Introduction 

Fading  prediction  methods  for  mobile  radio  channels  were  shown  to  enable  adaptive 
modulation  AM  [1-3]  in  narrowband,  Orthogonal  Frequency  Division  Modulation  (OFDM), 
adjacent  frequency  and  Multiple  Input  Multiple  Output  (MIMO)  fading  channels  in,  e.g.,  [4- 
8,10,25-27].  In  this  letter  and  [9,11,24],  we  explore  AM  aided  by  the  Long  range  prediction  (LRP) 
for  Slow  Frequency  Hopping  (SFH)  spread  spectrum  mobile  radio  systems  that  employ  coherent 
detection  [13,15].  We  propose  to  predict  the  channel  coefficients  in  the  next  hopping  frequency  of 
SFH  systems  based  on  a  number  of  past  fading  observations  from  previous  hopping  frequencies. 
Fading  prediction  is  challenging  in  this  case  since  past  observations  are  at  different  frequency  slots 
constrained  by  the  hopping  pattern.  Moreover,  we  investigate  joint  adaptive  transmission  that 
combines  frequency  diversity  and  AM  to  mitigate  the  effects  of  partial-band  interference  and  fading 
in  SFH  systems  with  coherent  detection. 

2.  System  Model  and  Long  Range  Prediction  for  SFH  Channels 

Consider  the  SFH  system  that  employs  coherent  detection  [15,16]  with  the  total  number  of 
frequencies  q,  the  hopping  rate  /„  and  the  frequency  separation  between  adjacent  carrier 
frequencies  A f.  In  this  letter,  we  employ  a  randomly  chosen  periodic  hopping  pattern  with  length 
N=q,  although  the  proposed  methods  also  apply  to  non-periodic  hopping  patterns.  Let  c((t)  be  the 
equivalent  lowpass  complex  sample  of  the  flat  fading  channel  at  time  t  and  frequency/  where /is 
the  carrier  frequency  (slot)  occupied  at  time  t  [13].  The  channel  coefficient  c(f,t)  is  closely 
approximated  by  a  zero  mean  complex  Gaussian  random  process  with  Rayleigh  distributed 
amplitude  and  uniformly  distributed  phase  [12],  and  we  assume  E\c(ft)\  =1.  The  spaced-time 

spaced-frequency  correlation  function  with  the  time  difference  x  and  the  frequency  separation  A f  is 
defined  as  [7,12,13,22]: 

R(Af,T)=E[c(f,t)c*(f+Af,t+x)]=  Rt(x)Rj(Af),  (1) 

where  the  factors  in  the  last  expression  are  the  time  and  frequency  correlation  functions. 

Figure  1  illustrates  the  adaptive  transmission  aided  by  the  LRP  for  this  FH  system.  Past  reliable 
observations  from  all  frequencies  are  fed  back  from  the  receiver  to  the  transmitter.  The  transmitter 
employs  the  LRP  to  predict  future  Channel  State  Information  (CSI),  and  adapts  the  transmission 
parameters  to  the  channel  variation.  We  employ  the  Minimum  Mean  Square  Error  (MMSE)  linear 
prediction  (LP)  algorithm.  Assume  the  channel  coefficients  c(f,t)  are  sampled  at  the  rate/=//7],  and 
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for  an  integer  n,  define  c(f(ji),ri)=c(J{nTs),nTs).  The  prediction  c(j{n+x),n+x)  (t  is  a  positive  integer) 

of  the  future  channel  coefficient  c(/(«+t),«+t)  based  on  p  past  observations  c(/(n), c(f(n- 
p+l),n-p+l)  is  fonned  as  (see  Figure  lb) 

c(J[n+T),n+T)=Xdj(n)c{f[n-j),n-j)  (2) 

j=o 

where  dj(n)  are  the  filter  coefficients  at  time  n,  and  % Ts  is  the  prediction  range.  Note  that  the 

sampling  rate  in  (2)  is  much  slower  than  the  symbol  rate,  but  faster  than  the  hopping  rate  fiu 

The  objective  is  to  find  the  LP  coefficients  that  minimize  the  MSE,  defined  as 

E [ | e(« ) 1 2]= F [ | c(f(n+i ),n+x )-c(f(n+i ),n+x)\  ’] .  Because  the  hopping  pattern  is  a  random  frequency 

sequence,  a  single  prediction  filter  does  not  exist,  and  the  LP  coefficients  need  to  be  re-computed  at 
the  sampling  rate.  The  optimal  LP  filter  used  at  sampling  time  n  is  given  by  [13]  d{n)=R{ny1  r(n), 

where  d(n)=[do(n).  ,.dp.i(n)]r,  R(n)  is  the  pxp  autocorrelation  matrix  of  the  observations  at  time  n 
with  Rjj(n)=E{c(J(n-i),n-i)c*(f(n-j),n-j)},  and  r(n)  is  the  cross-correlation  vector  of  the  observations 
and  the  prediction  at  time  n  given  by  rj(n)=E{c(/[n+x),n+x)c*(J[n-j),n-j)},  i,j=0,l,...,p-l,  for  given 

t.  Both  R(n)  and  r(n)  are  detennined  by  the  correlation  function  (1).  The  effect  of  additive  noise  on 
the  observations  can  be  incorporated  in  (2)  [4],  In  this  letter,  we  assume  that  the  noise  in  the 
observation  samples  is  negligible.  In  practice,  noise  reduction  techniques  can  be  employed  to 
improve  the  accuracy  of  the  prediction  [4,26,27].  For  realistic  SFH  systems,  the  prediction  MSE 
loss  is  dominated  by  observations  constrained  to  the  hopping  pattern  [9,1 1,24],  and  the  degradation 
due  to  additive  noise  is  relatively  small  in  this  MSE  region  [27]. 

For  realistic  mobile  radio  channels,  the  correlation  functions  /?,(t)  and  R/(Af)  in  (1)  must  be 

estimated  and  updated  when  new  observations  become  available.  We  employ  pilot  symbols  for 
estimation  [11].  The  rate  of  update  of  the  correlation  function  estimates  and  the  computational  load 
of  the  estimation  is  low  for  realistic  mobile  radio  channels  [7,11],  On  the  other  hand,  the  optimal 
MMSE  channel  prediction  method  (2)  is  complex  (on  the  order  of  p  multiplications  [28]),  because 
it  requires  inversion  of  a  large  matrix  at  the  sampling  rate.  To  reduce  the  complexity  to  the  order  of 
p“,  while  maintaining  the  same  performance,  we  employ  a  recursive  procedure  for  updating  this 
inverse  as  described  in  [9,11].  While  it  is  possible  to  reduce  complexity  further  by  employing,  e.g., 
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the  simplified  LRP  method  in  [9],  we  have  observed  that  suboptimal  prediction  methods  greatly 
degrade  performance  of  adaptive  FH  systems. 

3.  Adaptive  Modulation  Aided  by  LRP 

LRP  is  employed  to  enable  AM  for  each  upcoming  dwell  interval.  We  employ  adaptive  discrete 
power  discrete  rate  MQAM  with  M=2,  22(1_1),  i=2,3,4  [1],  As  in  [3,24,26],  reliable  performance  is 
maintained  by  incorporating  the  accuracy  of  the  predicted  CSI  into  the  AM  design.  The  symbol  rate 
is  20Ksps  (symbols  per  second),  and  the  target  Bit  Error  Rate  (BERtg)  is  10'  .  The  SNR  is  defined 
as  the  ratio  of  the  average  transmitted  symbol  power  Es  to  the  complex  white  noise  power  spectral 
density  N0  [1].  While  perfonnance  with  and  without  prediction  is  usually  compared  in  other 
adaptive  systems  (e.g  narrowband  or  OFDM  [4,8]),  this  comparison  is  meaningless  in  adaptive  SFH 
systems.  The  delay  associated  with  the  feedback  and  other  system  constraints  is  comparable  with 
the  dwell  interval  duration,  so  the  channel  estimates  obtained  during  current  dwell  interval  cannot 
be  employed.  A  single  outdated  estimate  at  the  previous  hopping  frequency  is  not  helpful  for 
enabling  AM.  Thus,  the  only  practical  alternative  to  using  fading  prediction  for  SFH  is  to  resort  to 
non-adaptive  modulation.  Therefore,  we  compare  the  Bit  Per  Symbol  (BPS)  of  AM  enabled  by  the 
LRP  with  that  of  non-adaptive  MQAM. 

In  the  numerical  results,  we  employ  the  standard  Jakes  fading  model,  a  typical  slow  hopping 
rate  fh=500  hops/second,  the  number  of  frequencies  q=32,  and  the  feedback  delay  of  at  least  1ms. 
To  enable  AM,  several  samples  are  predicted  for  each  dwell  interval,  and  the  average  prediction 
range  is  2ms.  The  maximum  Doppler  shift  is  50Hz  (equivalently,  the  prediction  of  0.1  wavelengths 
ahead  is  illustrated).  In  LRP,  we  employ  the  near-optimal  sampling  rate  2  kHz  and  p=50  in  (2)  [11]. 

In  Figure  2,  the  spectral  efficiency  (BPS)  of  AM  vs.  average  SNR  is  shown.  We  observe  that 
significant  gain  can  be  achieved  relative  to  non-adaptive  modulation  (Binary  and  Quaternary  Phase 

Shift  Keying  (BPSK,  QPSK)).  The  gain  depends  on  the  normalized  frequency  separation  Afa  given 
by  the  product  of  the  frequency  separation  between  two  adjacent  hopping  frequencies  and  the  rms 
delay  spread  a  of  the  fading  channel  [12-14].  In  [9,1 1,24],  we  have  demonstrated  that  the  prediction 

MMSE  increases  as  Afa  grows  since  the  observations  and  the  prediction  become  less  correlated. 
This  dependency  affects  the  bit  rate  of  AM  as  shown  in  Fig.  2.  For  small  A/a=0.01,  the  BPS  with 
prediction  approaches  that  with  perfect  CSI  [8].  While  the  BPS  diminishes  as  A/a  increases,  even 
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for  large  A/g=0.1,  the  gain  relative  to  the  non-adaptive  modulation  is  about  3dB,  or  1  BPS.  A 
physical  model  proposed  in  [4-6]  was  used  in  [11,24]  to  investigate  the  performance  of  AM  for 
SFH  systems  in  realistic  fading  channels.  While  the  BPS  is  lower  for  the  physical  model  than  for 
the  Jakes  model  due  to  the  time -variant  correlation  function  (1),  it  has  been  demonstrated  that 
significant  improvement  is  still  achieved  relative  to  non-adaptive  modulation. 

We  observe  that  the  FH  system  benefits  from  adaptive  transmission  primarily  when  A fa  does 

not  significantly  exceed  0.1.  Suppose  a  is  lps,  representative  of  suburban  areas  [12].  Then  a  SFH 
system  would  benefit  from  adaptive  transmission  when  the  frequency  separation  is  as  large  as  100 
KHz  (A/g«0.1).  In  realistic  SFH  systems  [13-16],  the  symbol  rate  is  on  the  order  of  tens  Ksps,  and, 
thus,  adaptive  transmission  aided  by  the  proposed  channel  prediction  method  is  feasible  for  these 
systems.  More  generally,  we  have  found  that  AM  is  feasible  in  SFH  systems  when  fdm<100Hz,  and 

the  total  normalized  bandwidth  (TNB)  qAfo  is  on  the  order  of  3  or  lower,  or,  equivalently,  the  total 

bandwidth  qAf  does  not  exceed  approximately  15  times  the  coherence  bandwidth  Bc«  1/5g  [14].  As 

the  TNB  grows,  the  spectral  efficiency  of  AM  saturates  and  approaches  that  of  non-adaptive 
modulation.  On  the  other  hand,  frequency  diversity  is  usually  exploited  in  FH  communications 
[13],  and  its  benefit  increases  as  the  TNB  grows.  Thus,  adaptive  transmission  and  diversity 
combining  compliment  each  other  over  the  practical  range  of  frequency  correlations  in  SFH 
systems. 

Our  results  demonstrate  that  fading  prediction  is  less  accurate  for  SFH  systems  than  for 
narrowband  transmission  [4,26],  OFDM  [8],  direct  sequence  CDMA  [4]  and  even  when  the 
observations  are  at  an  adjacent  frequency  [7].  This  loss  is  due  to  the  fact  that  the  observations  are 
constrained  by  the  hopping  pattern  in  LRP  for  SFH  systems,  and,  thus,  are  widely  distributed  in 
frequency.  This  constraint  degrades  prediction  accuracy,  and  hampers  utilization  of  fast  and 
efficient  adaptive  tracking  techniques.  However,  we  note  that  fading  prediction  is  critical  in  SFH 
applications,  since  adaptive  transmission  would  not  be  possible  without  prediction  in  FH  systems. 

4.  Adaptive  SFH  Systems  with  Partial-band  Interference 

We  focus  on  the  Partial-Band  Interference  (PBI)  that  is  not  due  to  a  hostile  jammer  [13,17,18], 
It  is  usually  slowly  varying  and  modeled  as  narrow-band  additive  Gaussian  noise  with  the  average 

power  spectral  density  Ni  with  that  occupies  a  small  fraction  8  of  the  total  bandwidth  of  the  FH 
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system.  We  use  adaptive  frequency  diversity  to  jointly  mitigate  the  effect  of  PBI  and  fading.  In  the 
proposed  diversity  FH  method,  the  same  information  is  transmitted  on  several  carrier  frequencies 
chosen  according  to  a  hopping  pattern,  and  the  outputs  of  different  diversity  branches  are  combined 

at  the  receiver.  For  simplicity,  in  this  letter  we  employ  only  two  frequencies  (diversity  branches)  / 1 

and  /  with  large  separation  NAf/2  and  negligible  correlation  assumed  assured  by  the  appropriate 
hopping  pattern  design  [11,20], 

We  assume  that  the  receiver  knows  perfectly  where  the  PBI  is  present  [11,13,19].  At  the 
transmitter,  the  uncertainty  of  PBI  presence  for  two  upcoming  frequencies  f1  and  j  is  modeled  as 

k 

follows.  Define  the  indicator  function  for  the  presence  of  PBI  at  the  upcoming  frequency/  ,  k=l,2, 

k 

as  /  =1  if  the  interference  is  present  at/  ,  and  0  otherwise.  Given  the  reliability  factor  qe[0,/],  the 
probability  of  the  interference  at  the  transmitter  is  modeled  as 

Pk=  r\h+(Fr\)(l -  Ik).  (3) 

The  following  near-optimal  diversity  combining  technique  [11,13,21]  is  employed  at  the 
receiver.  When  there  is  no  interference  at  both  hopping  frequencies  /  and  /  ,  Maximal  Ratio 
Combining  (MRC)  is  used.  When  only  one  frequency  has  interference,  the  PBI-free  branch  is 
selected  using  Selective  Combining  (SC).  If  both  frequencies  are  interfered,  random  guess  is  used 
to  detect  data.  More  complex  sentient  FH  technique  (similar  to  the  selective  transmitter  diversity 
[10])  can  further  improve  performance.  In  this  method,  channel  coefficients  of  L  widely  spaced 
frequencies  (chosen  according  to  the  hopping  pattern)  are  predicted,  and  a  subset  of  r  frequencies 
with  the  largest  channel  gains  are  selected  in  the  transmission.  Due  to  lower  average  transmitted 
power  [11],  sentient  FH  also  results  in  lower  Multiple  Access  Interference  (MAI)  than  FH  given  the 
same  number  of  transmitted  frequencies  r. 

The  degrading  effect  of  PBI  is  compensated  for  by  the  benefit  of  frequency  diversity  on  the 
performance  of  LRP  for  these  diversity  systems,  and  the  prediction  accuracy  is  better  than  for 

interference-free  systems  without  diversity  [11],  In  AM  for  diversity  FH,  let  at  and  a  a ,  k=l,2,  be 

the  actual  and  predicted  fading  channel  amplitudes  at  the  upcoming  frequencies  /  and  /  , 
respectively.  The  average  BER  of  fixed  power  AM  when  M(z')-QAM  is  employed  by  the  transmitter 
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*  «>  «>  E  (  Ctl  ) 

BERW„  (Es,N0,ai,a2,pi,P2)  =  (l-pi)(l-p2)f  jBER,w/y  =  ~  N — —)p( a / 1 a , )p{a2\a2 )da  /da? 

0  0 

oo  ^7  00  Ej  C%2 

+ p2(I-pi)j BERiW/y=  SN‘  )p(ai\af)&ai  +  /?/(i-/?2)  jBERM//y=-^)/?(a?|a?)da?  +  0.5 p2p2,  (4) 

0  0 

where  the  BERM^  is  an  upper  bound  on  the  BER  of  M(i)-QAM  [1],  /ha.///)  is  the  conditional  pdf 
of  a*-  given  </  [3],  and  pk  is  given  by  (3),  k=l,2.  The  modulation  level  is  chosen  asM=max{M(z)  | 
BER M(i)*(Es,No,&i,a2,  P1P2)  ^  BERtg}.  In  AM  combined  with  sentient  FH,  for  r=2,  (4)  is  employed, 

assuming/7  and  j  are  the  two  frequencies  (among  L)  with  the  largest  prediction  gains. 

Simulations  are  used  to  demonstrate  the  performance  of  adaptive  SFH  for  typical  PBI  values 
[23].  In  Figure  3,  the  BPS  of  AM  is  illustrated  under  the  assumption  of  perfect  knowledge  of  PBI  at 

the  transmitter  (r)=l).  We  observe  that  the  BPS  degrades  as  5  increases,  and  sentient  FH 

outperforms  diversity  FH.  As  L  grows,  the  BPS  of  sentient  FH  saturates,  and  L= 4  has  near-optimal 
performance  [11].  We  also  show  that  the  BPS  of  adaptive  modulation  for  a  non-diversity  system 
(single  Rayleigh  fading  channel)  with  PBI  is  poor  [11].  (Note  that  when  this  method  is  extended  to 

r)  <1,  the  target  BER  cannot  be  satisfied,  implying  that  diversity  is  required  for  channels  with 
imperfect  knowledge  of  PBI  at  the  transmitter.)  Figure  4  shows  the  BPS  degradation  of  diversity 
FH  as  t]  decreases.  For  both  diversity  and  sentient  FH,  when  r\<0.95,  the  target  BER  cannot  be 

satisfied  with  the  adaptive  transmission  method  proposed  above,  and  additional  diversity  would  be 
required  to  maintain  reliable  performance. 

5.  Conclusion 

The  optimal  MMSE  long  range  prediction  algorithm  for  SFH  communications  with  coherent 
detection  was  introduced.  It  was  demonstrated  that  that  the  proposed  LRP  method  enables  adaptive 
modulation  for  SFH.  Moreover,  joint  adaptive  frequency  diversity  and  AM  was  investigated  for 
channels  with  PBI.  Numerical  and  simulation  results  demonstrate  that  significant  performance 
gains  can  be  achieved  relative  to  non-adaptive  modulation  for  realistic  FH  systems. 
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b.  The  LRP  algorithm 

Figure  1.  Adaptive  transmission  for  FH  channels  aided  by  long  range  prediction. 


Figure  2.  Spectral  efficiency  of  adaptive  modulation  using  Long  Range  Prediction, TFs=2ms, 


fdm=5(Mz,  q=32. 
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Figure  3.  Performance  of  adaptive  SFH  with  partial-band  interference,  Jakes  model,  prediction 
interval  x Ts=2m s ,  fdm=50 Hz,  A fo=0.05,  8=0.1,  p =1.0,  Es/Ni=0dB,  BERtg=10'3,  L=4,  r=2  for  sentient 


FH. 
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Figure  4.  Performance  of  adaptive  diversity  FFI  with  PBI,  Jakes  model,  8=0.1,  xTs=2ms, 
fdrn=50Hz,  Afa=0.05,  Es/Ni=0dB,  BERtg=10'3. 
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